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ABSTRACT 


Ue report on the development of a computer-assisted method 
for the determination of the angular distribution data for 
secondary particles produced in relativistic nucleus-nucleus 
collisions in emulsions. The method is applied to emulsion 
detectors that were placed in a constant, uniform magnetic 
field and exposed to beams of 60 and 200 GeV/nucleon 
16 q ions at the Super Proton Synchrotron (SPS) of the 
European Center for Nuclear Research (CERN). Linear 
regression analysis is used to determine the azimuthal and 
polar emission angles from measured track coordinate data. 
The software, written in BASIC, is designed to be machine 
independent, and adaptable to an automated system for 
acquiring the track coordinates. The fitting algorithm is 
deterministic, and takes into account the experimental 
uncertainty in the measured points. Further, a procedure 
for using the track data to estimate the linear momenta of 
the charged particles observed in the detectors is included. 
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I 


INTRODUCTION 


In this paper we report on the development of a data 
analysis method for the rapid determination of the azimuthal 
and polar emission angles of particles produced in nucleus- 
nucleus collisions observed in emulsion chambers exposed to 
relativistic ■*•^0 beams. The method makes use of the track 
coordinate (x,y,z) data that is presently obtained by visual 
inspection of the developed emulsion plates, using scanning 
microscopes. Although our initial application focusses on 
studying charged pions, the method is applicable to data for 
any emitted particle. The results of this work will be 
applied to the analysis of heavy ion cosmic ray interactions 
that are observed in emulsion chambers flown at high 
altitudes^. Events from these cosmic ray experiments are 
especially valuable since they often occur at energies that 
are substantially greater than those readily achievable with 
present-day particle accelerators. 


The angular distributions of secondary particles, generated 
in the collision of two nuclei, contain information on the 
dynamics of the nuclear interaction process. Events that 
are characterized by large numbers of secondary particles 
and large transverse momenta are likely candidates to 
exhibit new fundamental phenomena. One such phenomenon is a 
new state of matter, the quark-gluon plasma (QGP), that is 
expected to occur in relativistic collisions that involve 
unusually high energy densities 2 . Another example rests in 
the idea that, if the collisions are simple superpositions 
of proton-like collisions, the produced particles are 
expected to be emitted isotropically in the center-of-mass 
frame. In each of these cases, it is very important to 
examine and understand the angular distributions of 
particles produced in high energy nuclear interactions, 
specifically with respect to non-statistical structure that 
may contain signatures of new physics 3 ' 4 . 

Emulsion chambers are well established as a tool for 
observing nuclear interactions involving energetic charged 
projectiles. They have the advantages of being relatively 
durable and easy to prepare. They can be used to accurately 
measure the charge and energy of the primary projectile, in 
addition to the emission angles associated with fragments 
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and secondary particles produced for the highest energy 
nuclear interactions. However* being passive detectors* 
they require a lengthy and somewhat involved set of 
developing and scanning procedures in order to obtain the 
raw data needed for analyzing the events they record. Even 
after the emulsion plates are developed, considerable 
laboratory work is needed to obtain angular distribution 
data . 
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II. OBJECTIVES 


The primary objective of this project is to analyze 
secondary particle distribution data, recorded in emulsions 
from the EMU05 experiment 5 , for the existence of non- 
statistical structures. To accomplish this objective, it 
was necessary to develop appropriate computer software that 
could be used to find the azimuthal and polar emission 
angles from particle track coordinate data. The software 
includes error analysis, and it has been tested successfully 
with data for which the results are known. In particular, 
the angular distributions of charged pions, observed in the 
EMU05 experiment, are to be examined for deviations from 
isotropy in the center-of-mass frame. 
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III. 


EXPERIMENT DESCRIPTION 


For the EMU05 experiment, pulsed beams of with energies 

of 60 and 200 GeV/nucleon were provided by the Super Proton 
Synchrotron (SPS) at the European Center for Nuclear 
Research (CERN). The pulse duration was 2s with a total 
intensity of 3xl0 3 ions/cm 2 pulse. The integrated exposure 
given to a chamber was 10^ ions. The beam size was 
2 . 54cmx2 . 54cm (1 sq . in.) and each chamber was exposed to 
beam spills shifted laterally from each other by 1 cm. 
Proportional counter measurements at the chamber, located 30 
cm downstream from the beamline end, indicated the beam to 
be 98% pure. 

The chamber used in this work consisted of stacked emulsion 
plates separated by layers of lead, CR39 plastic and 
polystyrene. A sketch of the experimental arrangement, 
showing the approximate dimensions of the chamber, is 
provided in figure 1. The chamber was placed inside a 
uniform, 1.8 Tesla magnetic field. A cross sectional view 
of the detector configuration for which our analysis method 
was developed, is provided in figure 2. In this case, each 
emulsion plate had a 70ym base coated on both sides 
with 50pm of emulsion. The separation between the emulsion 
plates is not constant, but gradually increases in the 
direction of the beam. This facilitates the measurement of 
the track curvature, the identification of the charge of the 
emitted particle, and places an upper limit of 10 GeV on the 
energy of the secondary particles that can be analyzed. 

Also, lead plates are placed near the front of the detector 
where the density of emulsion plates is greater to increase 
the likelihood of collisions there. This feature also 
improves the accuracy with which the position of the 
collision vertex and the track angles can be determined. 
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IV 


DATA REDUCTION METHODS 


IV. a. Determination of Emission Angles 

The method employed to find the polar and azimuthal emission 
angles consists essentially of fitting the set of position 
coordinates, ( x i'yi/ z i)/ for a given track, to the equation 
of a straight line®. The situation is illustrated in figure 
3. The vector d points in the initial direction of motion 
of the emitted particle. Since the paths are curved, it is 
recognized from the outset that this approach can be used to 
obtain a good estimate of the initial direction of motion of 
the outgoing particle, at the point of collision. 
Consequently, only those points closest to the collision 
vertex are used in the calculation. 

First, a fit to the line y = a + bx is found using the set 
of points, (x i ,y i ), in the x-y plane. The azimuthal angle, 
0 , is then simply obtained from 


0 = tan 1 ( b ) 


( 1 ) 


where b is the slope of the line. The procedure is repeated 
for the set of points, (r^,z^), in the r-z plane where 


r . 

i 


/ 2 2 

u.r + (y t r 


( 2 ) 


and z = c + mr. The angle 0 is then obtained from 


0 


tan 


(m) 


(3) 


This procedure is performed for each track associated with 
the event. 

Values for b and m are obtained by the minimization of a 
chi-square quantity given by 
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X 2 (a,b) 


2 


( 4 ) 


N 

£ 

i = l 



y (x i ;a,b) 


) 


t 


where is the experimental uncertainty in the ith point. 
The resulting conditions. 



and 



(5) 


must be satisfied, and in doing so yield two equations in 
two unknowns that are readily solvable for the constants 
a(c) and b(m) . An estimate of the probable uncertainties in 
the constants can be obtained if the data are treated as 
independent with each contributing its own bit of 
uncertainty to the parameters. Consideration of the 
propagation of errors shows that the variance, a^, in the 
value of any function will be 


where f 


°f 


2 


N 

£ 

i = l 




2 


a (c ) , b (m) . 


( 6 ) 


If, however, the individual measurement errors of the points 
a^, are not known, then a more accurate estimate of the 
probable uncertainties in the parameters a(c) and b(m) can 
be obtained via the following procedure. Set = 1 in 
equations (4), (5), and (6), and multiply the values 
of cr^ by the additional factor. 


/ x 2 /(N-2) ' 

2 

where x is computed by (4). In essence, this latter 
procedure is equivalent to assuming that one obtains a good 
f it . 
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IV. b Calculations of the Linear Momenta 

As suggested in the introduction, it is important to 
identify those events in which a large amount of linear 
momentum of the incident projectile is transferred to the 
target nucleus. This may be done by careful examination of 
the linear momenta reaction products 7 . 

The radius of curvature of the path of a secondary particle 
is related directly to its linear momentum. To show this, 
consider the motion of a charged particle in a magnetic 
field. The magnetic force on the particle is given by 

F = q(vxB) , (7 ) 

where q is the charge on the particle, v is its velocity, 
and B is the magnetic field. In the present case, B is 
assumed to be uniform and oriented in the positive y- 
direction. Thus, the magnitude of the force can be written 


F = qvbsin ( 0 ' ) 


( 8 ) 


where 0' is the angle between v and B, and the direction of 
F is everywhere perpendicular to the plane formed by v and 
B. The curved motion is described in terms of a centripetal 
acceleration so that. 


qvBsin ( 0 ' ) = mv 2 /R , ( 9 ) 

where m is the mass of the particle and R is its radius of 
curvature. Since p = mv is the linear momentum of the 
particle, we have 

p = qBRsin 0 ' . (10) 

For convenience, equation (10) may be expressed as^ 

p(GeV/c ) = 0.29979 q B(T) R(cm) sin0' , (11) 

where q takes on the value ±1 for pions. 
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We can derive an estimate of R from the measured track 
coordinates using the scheme illustrated in figure 4. From 
the figure. 


= R sina ^ 


( 12 ) 


where L is the distance along the symmetry axis of the 
detector, in this case the z-direction, to the ith emulsion 
plate. Also, note that 


Ax ^ = R(1 - cosa. ) 


(13) 


The quantity Ax^ is the perpendicular distance from the beam 
direction. These last two equations can be combined to give 


R = Ax^/(.l - cos(sin 1 ( L ^/R ) ) ) . 


(14) 


Since we are interested in obtaining a solution to this 
last, non-linear equation for R in terms of Ax^,_£his may 
best be done by approximating the cos(x) and sm (x) 
functions by their series forms, i.e., 


cos ( x ) 



+ x < °» 


. -1 
sm 


(x) 


X 



1 • 3 *x 


5 


2 • 4 • 5 


+ 


x < 1 . 


Using only the first order terms, we obtain 

R = (L.) 2 /2Ax.. 


Clearly, this approximation is best suited for measurements 
involving the coordinates of the first few emulsion plates 
nearest the interaction vertex, and for reaction products 
with large p values. 


XXIV-8 


V. SOFTWARE DEVELOPMENT 


A computer code that makes use of the analysis methods 
described in the previous section was written for the 
Commodore AMIGA 9 computer. The code is written in BASIC and 
is designed to be machine independent. It is expected that 
the code will be executed under the BASIC interpreter 
supplied with the computer. For input, the program requires 
files that contain the track coordinate (x,y,z) data that 
have been obtained for each event by scanning the developed 
emulsion plates. At present, the program returns the 
corresponding angles, 6 and 0 for each track, and it also 
has a provision for estimating the linear momenta of the 
emitted particles, based upon the track radius of curvature 
and the magnetic field. Early tests, employing idealized 
track data, as well as actual track data from a few plates, 
indicate that the code is operating correctly. A current 
source listing, to be regarded as preliminary, is provided 
in Appendix A along with a logic diagram for the code. A 
detailed description will appear elsewhere, after 
finalization of the software. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


Initial development work on the computer software for 
determining the emission angles and estimating the linear 
momenta of particles emitted in nucleus-nucleus collisions 
observed in emulsions has been completed. The software has 
been tested successfully for correct operation using 
idealized track data and partial data from the EMU05 
experiment. Further testing of the code with complete track 
data for EMU05 events is recommended to confirm the accuracy 
of the calculations. 

Additional heavy ion experiments involving emulsion chambers 
of the 5A2 design are planned for the SPS accelerator. The 
first will employ a 32g beam and is scheduled for September, 
1987. Another will use a 208pb beam that is anticipated 
being available during the Fall of 1989. Also, the High 
Energy Astrophysics Branch of SSL has been involved over the 
past 10 years in a collaborative research program, the 
Japanese American Collaborative Emulsion Experiment (JACEE), 
the purpose of which is to study charge particle cosmic ray 
interactions in emulsion chambers flown at high altitude. 

To date, seven balloon flights have been conducted and data 
analysis has been completed for five of these. In view of 
the large amount of data anticipated to be available from 
these two efforts, it is recommended that an automated 
system of coordinate data recording be incorporated with the 
code development work presently underway in order to reduce 
the time between plate scanning and final analysis of the 
angular distributions. Such a system will be especially 
valuable for analyzing events having high multiplicities. 

The present method of calculating the emission angles will 
work best when data are available for a few closely spaced 
plates near the interaction vertex of the event. It is 
therefore of interest to explore alternative means of 
fitting the track data that make use of functions that 
better represent the curved path. An initial approach would 
include using higher order polynomial function 
approximations to the path, and finding the tangent to the 
curve at the interaction vertex. 
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10 REM****PR0GRAM TO CALCULATE THE AZIMUTHAL AND POLAR EMISSION*** 

20 REM ANGLES FROM THE EMU05 EXPERIMENT DATA.***** 

30 REM***************EMU05B . BAS************************** 

40 DIM A<4> ,X<50) ,Y(50) ,Z<50) ,SIG(50) ,R<50> 

50 XMAX=80000! :VMAX*80000! :ZMAX-47o80! 

60 ANGFAC* < 1801/3.14159) :PFACfc*l 
SO CONST - .29979: BFIELD * 1.8 
90 QP » 1 ! : ON - -1 ! 

100 REM*****READ IN THE TRACK DATA.************** 

110 REM*****I 1NITS SHOULD BE MICRONS********************* 

200 INPUT *'F ILENAME=» " , FILENAME 
210 OPEN FILENAME FOR INPUT AS 1 
212 LPRINT "DATA FILE - ", FILENAMS 
215 REM ******** 

217 LPRINT 

220 REM********** 

225 INPUT #1. NR AY!: 

227 LPRINT " NRAY- " ,NRAY7. 

228 REM****** 

230 PPR - 0»: PPL - 0 ! : PPRTOT-O ! :PPLT0T»0 ! 

240 fcr j:; - l TO NRAYN 

2^5 LPRINT " INPUT TRACK DATA FOLLOWS " 

246 LPRINT 
250 NPTSt: * 0 
260 FOR K% =* 1 TO 40 
270 NPTSt: =* NPTSt: * 1 
275 REM******* 

280 INPUT #1, X(K!4> ,Y<K7.> ,Z(K7.) ,SIG(K!4> 

205 R ( K3C ) - SQR ( X ( K% ) *X ( K% ) + Y(K%)*Y(K%) ) 

290 LPRINT USING "########.##"; X (KX> , Y ( ) ,Z ( K % ) , R ( KX ) , SIG < K%> 

300 IF X(K%) - -1 ! THEN NPTSX - NPTSX - 1: GOTO 320 
305 REM END INNER LOOP 
310 NEXT K7. 

320 REM*****PERFQRM FIT TO A STRAIGHT LINE AND****** 

330 REM ***BASEU ON THE FITTED DATA , FIND THE EMISSION ANGLES.***** 

332 LPRINT " EMISSION ANGLES FOLLOW » 

335 GOSUB 1000 

340 REM******************* 

345 GOTO 510 

350 REM*****BASED ON ITS ESTIMATED RADIUS OF CURVATURE, DETERMINE 

351 REM THE LINEAR MOMENTUM OF THE TRACK.********************* 

360 FOR IR7* * 1 TO NPTSt: 

370 R I - Z< IRK) *Z< IR7.) 

300 RI - RI/C2!*X< IR%)> 

390 R * R + RI 
400 NEXT IR7. 

410 FLN=NPTSt: : AVR = R/FLN: AVR=AVR/ 10000 ! 

420 PMOM » CONST *QP*BF I ELD*AVR 
430 PPR = PMQM*SIN< THETA) 

440 FPL = PMOM*COS( THETA) 

450 F PR TOT = PPRTOT + PPR 
460 PPLTOT =* PPLTOT + PPL 
470 REM*** ******* 

400 LPRINT USING "######### .##" ; PMOM , PPR , PPRTOT , PPL , PPLTOT 
490 REM**** 

500 LPRINT 

510 REM ****GET DATA FOR THE NEXT TRACK, OR 
530 REM END OUTER LOOP.********************* 

540 NEXT J% 

545 REM********************* 

550 CLOSE #1 

560 REM******************** 

570 END 
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ORIGINAL PAGE IS 
OF POOR QUALITY 



ORIGINAL PAGE IS 
5R0 REM OF POOR QUALITY 

1000 REM ***CUBPROGRAM TO PERFORM A BEST FIT TO STRAIGHT LINES; 

1005 REM IN THE X-Y AND R-2 PLANES. THESE FITS WILL BE********* 

1010 REM** TO DETERMINE THE AZIMUTHAL AND POLAR EMISSION ANGLES **-***"* 
o EM N0TEs ** Y ** B * x + A IS THE FORM OF THE STRAIGHT LINE . **■■***-* 
1015 REM*******F I RST CONSTRUCT THE SUMS OVER THE DATA POINTS 
10t-0 REM NEEDED FOR THE CALCULATION OF THE CONSTANTS. ****** 

1030 SX-0:SY-0:S=0:SXX»0:SXY-0 

1045 REM****USE ONLY A FRACTION OF THE AVAILABLE POINTS.***** 

1050 NPTSHS - NPTSS/PFACS 
1060 FOR J2S - 1 TO NPTSH*/. 

1070 SIGSQR - SIG( J27.)*SIG< J2S) 

1070 SXY » SXY + XC J2S)*Y< J2S) /SIGSQR 

1100 SXX - SXX + X« J2S)*X< J2S) /SIGSQR 

1120 SX - SX + X<J2S) /SIGSQR 

1130 SY - SY + Y(J2S> /SIGSQR 

1150 S - S * 1 ! /SIGSQR 

1160 NEXT J27. 


1 170 REM************************* 

1180 REM* ***THE NEXT STEP IS TO CALCULATE THE VALUES 

1190 REM OF THE PARAMETERS A AND B FOR THE AZIMUTHAL ANGLE ***** 

1FOO DELTA - S*SXX - SX*SX 

1210 APHI * ( SXX*SY - SX*SXY ) /DELTA 

It -20 BFHI - (S*SXY - SX»SY> /DELTA 

1230 REM *****FIND THE UNCERTAINTY IN THE A AND B.**** 

1 2*0 SIGMAA = SXX /DELT A ; S IGMAA = SQR(SIGMAA) 

1250 SIGMAS = S/DELTA; SIGMAS = SQR(SIGMAB) 

1760 REM**** * * * * * * ********* ********* * *** 

1770 REM*****CALCULATE A CHISQUARE VALUE FOR THE FIT.***** 

1730 REM AND RE-ESTIMATE THE UNCERTAINTIES IN A AND 3 *** 

1790 SRE3*0! 

1800 FOR IS - 1 TO NPTSH*/. 

1810 YP » BPHI *X (IS) + APHI 
1820 RESI = ( Y ( IS) - YP >/SIG(IS) 

1830 RESI * RESI*RESI 

1840 SRES = SRES + RESI 

1850 NEXT IS 

1060 CHISQR = SRES 

1870 EFACT - CHISQR/ (NPTSHS - 2) 

1880 EFACT - SQR( EFACT) 

1890 REM*****REEST I MATE THE UNCERTAINTY IN THE FITTED CONSTANTS ***** 
1900 SIGMAA = SIGMAA*EFACT 
1910 SIGMAS » SIGMAB*EFACT 

1912 REM*******NEXT, CALCULATE THE ANGLE PHI . *********** 

1915 PHI * ATNtBPHI ) :PHI = ANGFAC*FHI 

1917 LPRINT » A SIGMAA B SIGMAB PHI" 

1V20 REM********************************** 

1924 IF < X C 1 ) <0! AND Y(l) < 0!) THEN PHI ■ PHI + ISO'sGOTO 1°29 

1925 IF X(l> < 0! THEN PHI = PHI + 130! s GOTO 1929 

1926 IF Y < 1 ) < 0! THEN PHI = PHI + 360 ! s GOTO 1929 

1929 REM ************************ 

1930 LPRINT USING "######.##"; APHI , SIGMAA , BPHI .SIGMAB .PHI 
1940 LPRINT 


2000 REM* ***NOW DO THE SAME FOR THE POLAR ANGLE****** 

2005 REM * * ***X — > R AND Y — > Z. ********** 

2010 SR=0 : SZ=0 ; SRR=0 : SRZ=0 : S=0 
2005 FOR J1S - 1 TO NPTSHS 
2002 SIGSQR = SIG( J1S)*SIG( J1S) 

2010 SR = SR + R(J1S) /SIGSQR 

2015 SZ = SZ + Z( J1S)/SIGSQR 

2020 SRR = SRR + R « J 1 S ) *R ( J1 S ) /SIGSQR 

2030 SRZ = SRZ + R ( J1 S )*2<J1S > /SIGSQR 

2035 S - 5+ 1! /SIGSQR 

2040 NEXT J1S 

2042 REM * * ** *F I ND THE FITTED CONSTANTS, A AND B FOR THE DETERMINATION OF 

2043 REM THE ANGLES.***** 
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2043 DELTA - S*SRR - SR*SR 

2047 A THETA - ( SRR*SZ - SR*SRZ > /DELTA 

2050 B THETA =» <S*SRZ - SR*SZ) /DELTA 

2052 REM****N0W FIND THE UNCERTAINTIES IN THE FITTED CONSTANTS.******** 

2055 SRES = 0! 

2056 SIGMAA * SDR < SRR/DELTA > : SIGMAB - SQR(S/DELTA) 

2057 FOR 1% * 1 TO NPTSH7. 

2060 ZP = BTHETA*R( I’/.) + ATHETA 
2062 RES I = <Z(IV.) - ZP)/SIG(i:i) 

2065 RES I “ RESI*RESI 

2067 GRES = SRES + RES I 

2068 NEXT VA 

2070 CHISQR =* SRES 

2072 EFACT =* CH I SDR/ ( NPTSH*/. - 2) 

2075 EFACT = SQR( EFACT) 

2077 SIGMAA = S I GMAA *EF ACT 
2000 SIGMAB » SIGMAB*EFACT 

2081 REM ****** *NEXT, FIND THE ANGLE THETA.******************* 

2082 THETA =» ATN < 1 ! /BTHETA ) s THETA - ANGF AC *THET A 

2083 IF Z«l) O* THEN THETA - THETA *180! . 

2084 LPRINT " A SIGMAA B SIGMAB THETA" 

2085 LPRINT USING "######.##" JATHETA. S I GMAA » BTHETA » S IGMAB * THETA 

2086 LPRINT 

2090 REM * ******** ********************** 

2270 REM ************** 

2230 RETURN 
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